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depends strongly on the details of its thermal model, and the simplest theory for TES 
response assumes a single heat capacity connected to the heat bath. Here, analytical 
^ results are derived and discussed for the complex impedance, the responsivity and the 

1 — 1 noise of a transition-edge sensor, when the thermal model is not simple but consists 

of either two or three connected heat capacities. The implications of the differences 



of the models are discussed, as well. 
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I. INTRODUCTION 



Superconducting transition-edge sensors (TES) have become very popular as practical 

radiation detectors because of their high sensitivity and energy resolution, maturity of the 

fabrication technology required allowing wafer-scale processing, and versatility, so that they 

can be used for detection of both particles and photons from sub-millimeter frequencies 

to gamma-rays^. The basic operational theory of a voltage-biased TES in electrothermal 

feedback was described long ago'= however, in the original form the simplest possible thermal 

circuit was assumed, namely, that the detector could be described as a single lumped heat 

capacity, connected to the heat bath by a single thermal conductance. Sometimes this 

thermal model is fairly adequate in describing the detector response^, but some detector 

designs have been experimentally shown to behave in a more complex manner^]. Therefore, 

theoretical modeling has been advanced in recent years to include more complex thermal 
circuitsPraSl. 

One way to approach the problem is to generalize the problem fully to any number of heat 
capacity blocks and thermal conductances, and solve the obtained (large) linearized system 
of equations numericall y 11 * 16 !. This approach has the advantage that it is straightforward 
to move from simple models to more complex models within the same formalism, and for 
the most complex models this may be the only approach available. Nevertheless, for models 
of intermediate complexity, such as systems of two or three connected heat capacities, it 
is possible to calculate the detector properties analytically, as well. Analytical solutions 
have the advantage that they are easier to work with, can be used further for fitting of 
experimental data, and a lot of results can be calculated quickly. One can quickly compute, 
for example, how the detector response changes as a function of a certain parameter of the 
thermal model. 

Here, we present analytical solutions for the complex impedance, responsivity, and all 
the unavoidable components of the current noise spectral density of a TES detector, for 
all possible two-block thermal models and two simplest three-block models, in a compact 
formulation. We take advantage of the fact that we also derive the analytical equation for 
the complex impedance of the detector^CS Q f each model. Measurement of the complex 
impedance of a TES detector has been shown in recent years^^H to be a very valuable tool 
for characterising TES detectors; in particular for the discussion here, it gives information 
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FIG. 1. All possible two-block models studied, (a) Hanging model, (b) intermediate model, (c) 
parallel model. 



about the thermal model. There is some previously publishe d 4 * 10 * 12 * 15 '* and unpublished^ an- 
alytical work on certain two- and three-block models of bolometers. Those results describe 
some limits and special cases of this work. In RefsP^, some two-block models were dis- 
cussed, and an approximation was made about the steady state temperature of the blocks. 
In R e f s .112112 equations for the impedance, responsivity and the noise equivalent power were 
given for one two-block and three three-block models, but the non-Ohmic behavior of the 
bolometei^ZI was no t fully discussed. In RefP^, two simple three-block models were dis- 
cussed, but no simple analytical equations for the noise were given. The goal of the present 
work is thus to give an extensive set of equations for the responsivity, current noise and 
impedance in a compact and usable form, which can easily be used to analyze real noise 
and impedance data of TES detectors^. We also give many example plots to show how 
various thermal parameters affect the detector properties. Detailed discussion on the noise 
equivalent power, energy resolution and other figures of merit are left for future publications. 



II. TWO-BLOCK MODELS 

We discuss here all three possible two-block cases, shown in Fig [I] One heat capacity 
block is always the TES film, where the Joule power is dissipated, whereas the second block 
describes an additional thermal body, which could represent for example the insulating 
membrane, the absorber etcP^*^. The formalism is kept general so that we do not need to 
decide on the physical picture a priori. 

The first case, shown in Fig. [I] (a), is where the extra heat capacity C\ is connected 
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to the TES heat capacity Ct es , but nowhere else. This type of model has sometimes been 
called the hanging or dangling heat capacity model. In the second case, shown in Fig. [I] (b), 
the extra heat capacity lies in between the TES and the heat bath. The main difference to 
the first model is that now the steady state temperature of the extra block is not equal to 
the TES temperature, unlike in the first case of Fig. [I] (a). We call this second model the 
intermediate model. Finally, in the last model in Fig. [I] (c) (numerically studied in Ref.^), 
the extra block is connected to both the TES and the heat bath as in case (b), but an 
additional parallel heat path exists between the TES and the bath. We call this the parallel 
model. 



A. Equations to be solved 

The starting point for all calculations are the differential equations that describe the 
heat balance in the thermal circuits and the electrical equation of the simplest Thevenin- 
equivalent input circuit of the TES, consisting here of an equivalent voltage V b , an equivalent 
resistance Rl, and an inductance L in series. For a simple one-block thermal circuit, there 
is only one thermal equation, so that one needs to solve a coupled set of only two differential 
equations, as reviewed in RefP. They are coupled because the resistance of the TES detector 
R(T,I) depends on temperature and current, and because the Joule heating power Pj in 
the TES depends on the electrical parameters. It is this coupling that produces the negative 
electrothermal feedback in the response of a voltage biased superconducting detector^. For 
a two-block thermal model, one additional differential equation is added to represent the 
extra thermal body. The exact form of the thermal equations naturally depends on which 
model in Fig. [I] one is analyzing. Here, we only show the equations for case (a), the hanging 
model, other cases follow analogously. The set of equations, excluding the noise terms for a 
moment, reads 

Ci^ = A(T? es -T?), 

L^ = V b -IR L -IR(T tes ,I), (1) 
at 

where T tes and T\ are the instantaneous temperatures of the TES and the extra heat 
capacity, respectively, and the heat bath is at a temperature T b . A and B are constants 
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describing the strength of the thermal links between the TES and Ci, and TES and heat 
bath, respectively, and n and m are temperature exponents, which depend on the physical 
nature of the thermal linkdP. 8P, on the other hand, is some time-dependent power input 
to the TES. Using these equations, the full time-dependent response (temperature and 
current) could, in principle, be solved. However, the power flows are non-linear functions 
of temperature, and the TES resistance is a non-linear function of both T and /, so that 
an analytical solution for the system is not tractable in the general (large signal) case. In 
this work, however, we are only interested in the impedance, noise and responsivity, all of 
which are traditionally defined only in the small-signal limit, where the non-linearities can 
be ignored. The linearization of Eqs. [T]is done by the usual Taylor expansions around steady 
state values^ (denoted by subscript 0): 

IR{T tes , I) « I R + aI {R /T tesfi )AT tes + (1 + (3)R AI, 

Pj « P + [2I R + 0{P o /I o )]AI + a(P /T tesfl )AT tes , 
A(T t n es - T") » A{T? esfi - T- ) + nAT^AT tes - nAT^AT,, 
B(T t Z - If) « B(7£ fl - Tn + mBT t ZjAT tes , (2) 

where AT tes = T tes -T tesfi , AI = I-I , AT X = Ti-T 1>0 , and a = (T tesfi /R )(dR/dT)\ Io and 
(3 = (I /R )(dR/dI)\T te3 o are the dimensionless (logarithmic) transition sensitivity param- 
eters. Substituting the above expansions into Eqs. [T] lead to the linearized set of equations 

C t es d ^]} es = (2I0R0 + P^r) AI + ( a— — g tes ,i - 9tes,b] AT tes + g^ATi + 5P, 

at \ IqJ \ T tesfi J 

C ^ d ^? L = 9tes,\ (AT ies - AT X ) , 
at 

L ^T = 6V b - AIR l ~ «^°AT tes - (1 + P)R AI, (3) 

at J-tesfi 

where we have defined the differential thermal conductances gtes,i — n AT™ e ~l and gt es ,b = 
mBTJ^~Q, and used the fact that for the hanging model Ti = T tes - From now on, we also 
simplify the notation as T = T tes . Also, after linearization, one can ignore the temperature 
dependence of the heat capacities and consider them to be constants [C tes = Ctes(T ), 
C\ = C\(Tq) ], as the corrections are in second order^. 

Eqs. [3] have been analyzed in time-domain in ReP^, with focus on pulse response and 
electrothermal stability analysis. Here, we proceed to focus on the frequency domain (also 
discussed in Ref.^for this model), which is natural for impedance and noise analysis. Fourier 



transforming Eqs. [3] and simplifying notation, one gets 

I R (2 + P) 1 1 ^ 

n 

vjLI u = V u - [R L + Ra(l + 0)} h ~ Ch{9 ^\ + 9tea ' b) T Uttea , (4) 

where we have denoted the Fourier amplitudes of all variables as X u , have defined T tes ,i = 
Ctes/gtes,i, and have defined two other important time constants and a dimensionless quantity 
Ch that appears in the place of the loop gain of the simple model as 

Ctes 

(gtes,i + gtes,b)(l —Ch)' 

Tl 



gtes,l 



Ch = ( ^ - (5) 

(<?tes,l + gtes,b)To 

Note that the definition of this "effective loop gain" Ch here is different from the standard 
one-block model one 11 /: = P a/ '(& es ,&T ), and therefore C H = g tes ,b/(gtes,i + gtes,b)£- In fact, 
if one calculates the true zero frequency, perfect voltage bias loop gain starting from the 
definitions®', one arrives at C even for this hanging model. This is intuitive, as the coupling 
to bath is the same. However, the new definition does have the benefit that for each different 
model (hanging, intermediate or parallel), the equations [4] have the same mathematical form, 
only the definition of the effective loop gain and the time constants change, as will be seen 
later. It is, of course, still possible to work with other, less compact notations, as was done 
for example in Refs J 13 1 15 ^. 



B. Complex impedance 

1. Hanging model 

The first quantity we want to derive from Eqs. [4] for the hanging model is the frequency 
dependent complex impedance of the TES, Z tes ,H- As the full circuit impedance Z circ is 
calculated by Z C i rc = V^/I^, we can subtract from it the known impedance of the circuit 
outside of the TES to define Z tes = Z circ — R L — iuL for the simplest bias circuit case, or 
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tes — 



arc 



— Zbi as in general. We can set = 0, and obtain 



Z tes , H = R (1+P) + - Ch —R (2+P) /\l + i UTl - 9 J^- _ *—\ , (6) 



with definitions of £#, tj and T\ given in Eqs. |5j This can be compared with the result for 
the simple one-block modeP 



R (i + P) + 



i?o(2 + /3) 



1 



(7) 



tes 



1 -£ 



1 + iwr/ 



where 77 = Ct es /[<7t eSi &(l — £)] for the one-block model. We see that due to the second heat 
capacity, a new frequency- dependent term with time constant T\ appears in the denominator. 
Thus, the effect of C\ is non-linear. The strength of the extra term is not only set by the 
value of Ci, but also by how large the thermal conductance gt es ,i is relative to gt es ,b- However, 
as the effective loop gain Ch also depends on the ratio of the two thermal conductances, it 
turns out that Z tes ,H deviates maximally in the complex plane from the simple model for 
some value a = gtes,i/(gtes,i + 9tes,b) < 1, depending on the value of loop gain. In the limits 
a —7- and a — > 1, Z teSt H approaches the simple one-block model result, with a heat capacity 
C tes and C tes + Ci, respectively. We show some example plots of the effects of C\ and gt es ,i 
in Fig. [2} We see that increasing C\ distorts Z tes ^ more from the simple-model half-circle 
[Fig. [2] (a)], but keeps the direction of the extra "bulge" constant in the complex plane. 
Changing gt es ,i, on the other hand, changes both the location and the size of the "bulge" 
feature. The loop gain also has a direct and quite complex effect, as shown in Fig. [2] (c). 
In the large loop gain limit, the two-block hanging model approaches the simple model, but 
for typical values of C the effect of the extra term is strong, and, interestingly, remains even 
in the limit C — > 0. For < C < 1, the effect of the two-block model is opposite compared 
to C > 1, in that the circular shape of the simple model is distorted inward. 

In Fig. [2] (d) we show two examples of a potential pitfall in analyzing impedance data. In 
the limits where gtes,i is either clearly smaller or clearly larger than gt es ,b, the shape of Z tes ,H 
is only distorted at the low-frequency or high-frequency end, respectively. In that case, an 
incomplete frequency range of the data could lead to a misinterpretation of a two-block 
model as a one-block model. For example, in Fig. [2] (d) we show a plot of a two-block model 
with a = 0.9, (5 = 1, C = 1.65, and C\jCtes = 10, which can be fitted with a one-block model 
up to a frequency ur tes ~ 2, but with misidentified = 1.5, £ = 1.76 and C = 12C tes - The 
same way, the two-block impedance with a = 0.1 (other parameters the same) can be fitted 
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FIG. 2. (Color online) Complex impedance of a two-block hanging model, with varying param- 
eters. We have used Rq = 0.1Q and f3 = 1 in all plots, and frequencies run typically between 
uT tes = 0.01. .100, where T tes = C tes /gtes,b- (a) Z teStH as a function of C\, with Ci/C tes = 
0.33,0.49,0.73,1.1,1.65,2.48,3.71,5.57,8.35,12.5, C = 1.65 and a = g te s,l/(9tes,i + 9tes,b) = 0.5. 
Increasing C\ corresponds to a larger deviation from the one-block model, shown as dashed line, 
(b) Z tes , H as a function of a ( g tes ,i), with a = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, C = 1.65 and 
Ci/Ctes = 10. Increasing a corresponds to the bulge feature moving from the left, low- frequency 
side to the right, high-frequency side. Dashed line is the one-block model limit (a — > 0). (c) 
Z teSiH as a function of £ = P a/(g tes>b T ), with C = 0.49, 0.73, 1.1, 1.65, 2.48, 3.71, 5.57, 8.35, 12.5, 
Ci/Ctes = 10 and a = 0.5. Dashed lines show the corresponding one-block results, (d) Comparison 
of two Z tes ,H impedance curves (points) (C\jCtes = 10, £ = 1.65 and either a = 0.1 or a = 0.9) 
with one-block model fits (red lines). For a = 0.9, the one-block fitted parameters are f3 = 1.5, 
C = 1.76 and C = 12Ct es , whereas for a = 0.1§ we get fitted values C = 1.485 and C = 0.9C^ es , 
with (3 = 1 kept fixed. 



with a one-block model with misidentified C = 1.485 and C = 0.9Ct es , down to a frequency 
uTtes ~ 0.05. The high-frequency distortions can get especially tricky, as measurements 
always have an electrical high-frequency cut-off somewhere. 



2. Intermediate model 

Turning to the second, intermediate model (Fig. [l](b)), we can follow similar derivation 
as for the hanging model. Omitting details, we find the same mathematical form for Z tes , 
but with different definitions of the effective loop gain and time constants: 



. ■ 9tes,l{Tl) 1 
1 + ZLOTj — — 

\gtes,i{Ti) + gi,b (1 - Cim) 1 + iun 



(8) 

where 

Ctes 



Tl 



gtes,i(To)(l - Cm) 
9u8,i(Ti) +gi,b 



Cim — J^vr- (^) 

gtes^y 1 oj J o 

Here, we must also explicitly define at what temperature gtes,i is evaluated, as the two 
blocks have different steady-state temperatures in this case (T for C tes and 7\ for Ci), and 
therefore there are two different values for g tes ,i that come into play: gtes,i(Ti) and gtes,i(T Q ). 
gi b is always evaluated at T 1; and is therefore not explicitly written in the above equations 



for simplification. gx^iTb) will only come into play when calculating noise (Sect. II D). 

Eq. [8] thus shows that the complex impedance of the intermediate model behaves quali- 
tatively exactly like the hanging model one, so that similar plots to Fig. [2] can be generated, 
but with different parameter values. It is therefore very hard, if not impossible, to distin- 
guish between the hanging and intermediate models based on fitting impedance data alone, 
as was pointed out in RefsP^^ already. 

In the intermediate model, it is clear that the true loop gain is really different from the 
simple model, unlike in the hanging case, as the extra block affects the DC response. One 
way to find the loop gain is to look at the zero frequency limit of Eq. [8j Doing that, one 
finds C = P a/(g e ffT ), where 

gtes,l[J-l) + 9l,b 



is an effective differential thermal conductance, which is a series combination of the two 
individual thermal conductances. Notice, though, how gt es ,i has to be evaluated at different 
temperatures in the nominator and in the denominator. The added complication for analysis 
is that if Ti is not directly measurable (typical case), then the DC I-V measurements alone 
cannot fix To or the parameters A, B,n and m, as there are too many unknowns in the general 
case. Only if one can set n = m based on physical assumptions (same thermal conduction 
mechanism), can one determine both T and g e ff from the I-V data. T\ naturally depends 
on the relative strength of the two thermal conductances such that if n = m, 



gtes,l(Ti 



n 



gi,b 



l/n 



9tes,x( T l) + #1,6 



Tl 



(11) 



3. Parallel model 

Finally, we discuss the third and final two-block model, where both blocks have conduction 
channels to the bath, the parallel model [Fig. [T] (c)]. Again, it is straightforward to derive 
equations analogous to Eqs. [4] and to solve for Z tes , yielding once again the same form 



Ztes,P = R (l + 



Ro(2 + P)/ 



1 + iOJTj 



gtes,l(T ) gtes,l(Tl 



[g tes ,i(T ) + g t es,b][gtes,i(Ti) + fl^&Kl _ £p) 1 + iuj n 



(12) 



with definitions 



a 



Tl 
Tl 

c P 



ten 



(g te s,i(T ) + gtes,b)(l - Cp) 

gtes,x(T x ) +g 1>b 
PqO. 

(gtes,l(To) + gtes,b)To 



(13) 



Comparing with the two previous models, the amplitude factor for the extra term is 
naturally more complex (depending on all gs), but 77 and Lp are the same as for the 
hanging model, whereas T\ is the same as for the intermediate model. Again, similar plots 
could be produced as for the two previous models, and thus distinguishing the parallel 
model from the hanging and intermediate models by fitting impedance data alone is quite 
hopeless. In the parallel model, the true loop gain is still C = P a/(g e ffT ), but now the 
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effective conductance is slightly more complex, as it is a combination of parallel and series 
components: 

gtes,i(To)gi,b n .s 

9eff = (T ,, n + 9tes,b- (14) 

gtes,l{J-l) + gi,b 

Naturally, similar analysis issues exist for this model as was discussed for the intermediate 
model. 

C. Small-signal responsivity 

Next, we turn from complex impedance to small-signal current responsivity, which de- 
scribes the frequency-dependent current response of the device to power input to the TES. 
Naturally, one can define responsivities for power inputs to the other heat capacity blocks 
as well, but here we concentrate only in the direct responsivity, as it is important for the 
noise analysis discussed in the next section. The responsivity for each variant of the two- 
block models can be calculated again from equations ffl by keeping the power term P w 
and by dropping the voltage bias modulation term V u . Then, defining the responsivity as 
sj(u) = Iuj/Plo, one gets a result that looks a lot like Z tes , as we start from almost the same 
equations. Therefore, for all models, we can write the responsivity as a function of Z tes and 
the total circuit impedance Z circ = Z tes + Z bias in a compact form as 

s ^ _ 1 Z tes - R (l + 0) 

R (2 + /3) ' v ' 



where for the simplest bias circuit we have Z circ = Z tes + Ri + iuL. Eq. 15 above has the 
advantage that if the TES impedance is measured, as is commonly done, the responsivity 
can be immediately calculated from it. 

The effect of the loop gain on the current responsivity is well known, with higher loop gain 
increasing the responsivity at frequencies below the effective thermal time constant (where 
the responsivity starts to roll-off), and moving the time constant to higher frequencies. This 
takes place regardless of the complexity of the thermal model. On the other hand, the effect 
of changing G\ and gt es ,i is n °t so self-evident. In Fig. [3j we plot how G\ and gt es ,i m the 
two-block thermal circuits affect the responsivity, using the hanging model as an example 
(other models behave in a similar manner). We plot only the magnitude |s/(o;)|, as the phase 
is irrelevant for noise considerations, which is the focus in this work. We see from the plots 
that for the usual case where the electrical cut-off frequency is above the thermal cut-off 
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frequency (low enough inductance L), the effect of the extra thermal block is to reduce the 
responsivity in the mid- frequency range starting from rj~ . Increasing C\ moves the first 
partial cut-off set by rf 1 to lower frequencies [Fig. [3] (a)], whereas increasing g tes ,i makes 
the intermediate frequency suppression stronger [Fig. [3] (b)]. 
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FIG. 3. Responsivity of a two-block hanging model, with varying parameters. In both plots, we 
have used R = 0.1Q, (3 = 1, and r ei = L/[R L + R (l + /3)] = 0.015r tes , where r tes = C te s/9tes,b- 
(a) \ai(u)\ as a function of d, with Ci/C tes = 0.33,0.49,0.73,1.1,1.65,2.48,3.71,5.57,8.35,12.5, 
C = 1.65 and a = gtes,i/(9tes,i + 9tes,b) = 0.9. Increasing C\ corresponds to a shift of the partial 
thermal cut-off to lower frequencies, and the development of a "knee" in the intermediate frequency 
range, (b) as a function of a (or gtes,i), with a = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, 

C = 1.65 and C\jCt es = 10. Increasing a corresponds to the decrease of responsivity in the 
intermediate frequency range. Dashed line shows the simple model effective thermal time constant. 

D. Thermodynamic noise 

In TES detectors, there are several possible mechanisms for noise. Some of them are 
directly associated with superconductivity such as phase-slip shot nois d 26 | 27 J , fluctuation su- 
perconductivity nois e 1 28 1 29 - 1 and vortex motion noise®'. In all these mechanisms, there are 
external parameters such as the material quality of the superconducting film or the mag- 
netic field that affect the noise, and thus these superconductivity related noise sources can, 
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in principle, be controlled to some extent. In addition, there is currently no conclusive ev- 
idence for a significant role of the superconducting noise sources being dominant in TES 
detectors. In this paper we therefore only concentrate on the two unavoidable noise sources 
in TES detectors: (i) thermodynamic fluctuations of heat within the thermal circuit (ther- 
mal fluctuation noise, TFN)P^, and (ii) electrical fluctuations of the resistive circuit elements 
(Johnson noise pZI2U22|. 

The thermal fluctuation noise and Johnson noise of a simple one-block bolometer were 
discussed long agcP^, and their effect on a TES detector in negative electrothermal feedback 
has been reviewed thoroughly in RefW TFN and Johnson noise properties of some two- 
and three-body models have also been discussed in R e f s jil£U ° l 13 l 15 l ; but none of those Refs. 
give very simple compact analytical formulas for noise, nor do they discuss all the models 
covered in this work. In this section, compact analytical results for the current noise are 



given for all three variations of the two block models (and in Sect. Ill for two three-block 
models), using the already derived results for the complex impedance and responsivity. 



1. Hanging model 

a. Thermal fluctuation noise For each thermal conductance in the model, there are 
power fluctuations on top of the steady state power that is conducted. These fluctuations 
can be taken into account in a straightforward manner by adding appropriate fluctuation 
power terms SPi into the constituting thermal equations in Eqs. [TJ which for the hanging 
model then become (ignoring external power input) 

Ctes^ = Pj~A (I* - Tfi - 5P teS)1 - B (T- - If) - 5P tes , b 

C 1 d ^ = A(T t n es -T?) + 5P tes , 1 . (16) 

Note how the sign has to be different for the 5P tes ,i terms in the two equations, because an 
increase in outflowing power from i — > j is equivalent with a decrease of the outflowing power 
from j — > i. Linearization and transformation into frequency domain, as before, thus lead to 
equations that resemble the equations for complex impedance and responsivity. After some 
algebra, one can write surprisingly simple results (disregarding correlations between the the 
two fluctuating sources) for the two TFN current noise spectral densities Si^u) = 
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associated with the two thermal links g^. 

i r |2 Ptes,b \Ztes — Ro{~\- + fi)\ 2 



L u\tes,b 



\Z c irc\ 2 Io Rq(2 + /3) 2 



\T I 2 -PfL.i \Z tes ~ -Rq(1 + /3)| 2 cu 2 rf 2 2 u^rf 

1 u ' ltes ' 1 |Z circ | 2 / 2 Rl{2 + I3Y l + c^r 2 ^.i^AWJI 1 + W 2 r 2' l A 'J 

where we have defined the power TFN spectral densities of the two thermal links as P^ es b and 
P^ es 1 , and the total TFN current noise spectral density is simply the sum of the above two 
components Sj : tfn(^) = \IJ\tes b + \Iu\tes i- The exac t form of TFN power spectral densities 
Pf depend on the details of the nature of the thermal links, and results have been derived 
for phonon transport in the fully ballistic or 1-D diffusive limits, as reviewed in Ref.^, or for 
electron-phonon or Kapitza thermal conductance^^. For all cases where the temperature 
difference between two adjoining heat capacity blocks is a step function (ballistic phonons, 
e-p interaction, Kapitza resistance), the TFN power spectral density is well approximated 
by 

i* = 2k B {g i , j {T i )Tf + //, J (7})7j). (18) 



Looking at Eqs. 17, we see that the first term is the usual TFN noise term ("phonon 
noise") that exists in the one-block circuit, as well. Its frequency dependence is now simply 
complicated by the more complex responsivity, and the \Iu)\ 2 es b noise spectra would look just 
like the responsivity plots in Fig. [3] On the other hand, the second, new noise term looks 
very different: It has significant weight only above the frequencies set by the time constant 
of the extra heat capacity t\. 

b. Johnson noise of the TES The Johnson noise of the TES has a slightly more complex 
influence on the thermal and electrical equations: the noise is electrical in origin, so there 
is a fluctuation term in the electrical part of the constituent equations. In addition, the 
dissipated bias power also fluctuates, affecting the thermal part of the circuit. One way to 
properly take both effects into account was discussed in Ref.^, by introducing the so called 
internal impedance matrix. We follow the same approach, which leads to a different Taylor 
expansion of the Joule power (only in terms of current), and to a expression in frequency 
domain Pj = Iq[(Ro — Rl) — iwL]I u , if the simplest bias circuit is assumed. That leads to 
a set of equations in frequency space 

_ (Ro-R L -iuL)I 1 1 rp 

^tes Tl{L — Lh) T tes ,l 
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iwT u i — — {T utes — T u i) 
n 

iu;LI„ = V U) tea ~ [RL + Ro{l + P)\ & ~ ^*!±±9*2£ T (1 9 ) 

with the notation the same as for Eqs. |4j and where V^jes is the Johnson voltage noise of the 
TES, which in equilibrium has the well known expression V u t ea = y/4k B T Q R , but which has 



been showrP^to depend on the parameter f3 in the first order as V Wt tes = y 4&BTb.Ro(l + 2/3). 
After algebra, one gets an expression for the Johnson current noise spectral density 



Si j(ou) = which is expressed in terms of the known TES complex impedance as 



\Iu\J 



2 _ u,te 



V, 2 + os Z tes + Rq |2 



7 

^circ 



(20) 



R 2 (2 + f3y 

where we have again used the definition of the circuit impedance Z circ = Z tes + R L + iuiL. 
This noise adds in quadrature to the TFN noise terms, as it is uncorrelated with them. 

c. Johnson noise of the Thevenin (shunt) resistor, or other external voltage noise The 
Johnson noise due to the shunt and parasitic resistances (Rl) does not influence the thermal 
circuit, and is simply 



V 2 

I j |2 _ v u,sh / 91 x 

\ 1 ^\sh ~ i y i 2 > \ ZL ) 

^ rirr. 



where Z circ = Z tes + Rl + icoL is the full circuit impedance and 

V* sh = Ak B T sh R L , (22) 

if all of Rl is at temperature T s h. If Rl consists of two parts, say shunt resistance R s h and 
parasitic resistance^ R par a, and those two parts are at different temperatures T s h and T para , 
then more accurately V^ sh = AksiTshRsh + T para R para ) . Any other external voltage noise 
would also contribute through Eq. [2Tj but with V UjS h given by the external voltage noise 
amplitude. 

d. Noise plots To gain some intuition on the noise, we plot here some examples how 
the different TES current noise components vary with the thermal parameters. In Fig. |4| we 
first show how noise depends on G\. When C\jC tes is greater than w 0.5, a clear bump in the 
noise spectrum develops in the mid- frequency range between 1/n and the effective thermal 
cut-off set by Ct es - The onset of the bump moves to lower frequencies with increasing C\, and 
it is, for the used ratio of thermal conductances a = gtes,i/{gtes,i + Qtes,b) = 0.5, dominated 
by the thermal fluctuation noise component of the hanging block |/ w |t eSi i, as the phonon 
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noise |/ w |tea,6 actually decreases in the mid-frequency range. The magnitude of the bump 
saturates to a value determined by the thermal conductance gt es ,i when C\ is large enough. 
Johnson noise also develops a minor step feature, and although small in Fig. |4| its relative 
strength compared to the TFN noise depends on the details, such as the value of loop gain 
and resistance. The shunt resistor noise actually develops a dip for the values of L and a 
used in Fig. |4j but for the parameter values used here, its magnitude is small. 




1E-3 0.01 0.1 1 10 100 1E-3 0.01 0.1 1 10 100 
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FIG. 4. (Color online) Current noise of a two-block hanging model, with varying C\/Ctes = 
0.33, 0.49, 0.73, 1.1, 1.65, 2.48, 3.71, 5.57, 8.35, 12.5 in linear-log (a) and log-log (b) scales. Increasing 
C\ corresponds to a shift of the onset of the bump (or dip in |io>|tes,6 and to lower frequencies, 

and an increase of the magnitude of |iu|t eS) i noise component. Other parameters: C = 1.65, 
a = gtes,i/{gtes,i+ gtesfi) = 0.5, R = O.m, p = 1, R L = O.OOlfi, J = 10//A, g t es,b = 1 nW/K, 
r e i = L/[Rl + Rq(1 + /?)] = 0.015r tes , r tes = Ct es /gtes,b- Dashed line shows the simple model 
effective thermal time constant. 

In Fig. [5j the effect of the value of the hanging thermal conductance gtes,i is studied, by 
plotting the noise spectra as a function of the varying relative strength a = gt es ,i/ (9tes,i + 
9tes,b)- The difference to the effect of C\ is that instead of developing a bump with constant 
magnitude, the bump magnitude evolves with a in a non-monotonous manner, having a 
maximum at around a ~ 0.6. In addition, the bump shifts up in frequency with increasing 
a (both low and high frequency sides). In other words, the effect of the hanging block on 
noise is largest for approximately equal magnitudes of the two thermal conductances, just 
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like for the complex impedance. Interestingly, the bump in total noise eventually develops 
into a step-down feature for the highest values of a here (a = 0.9), where in a certain 
region of frequencies the noise can be lower than the simple model noise. This does not 
imply, however, that the noise equivalent power NEP = \I u \tot/\si{w)\ is lower, because the 
responsivity is reduced even more than noise in that frequency region (See Fig. [3]). Looking 
at the different components of the noise, one sees that the non-monotonous behavior in total 
noise is caused by the hanging block TFN noise |/ u |te S) i, whereas the phonon noise \I u \t es ,b 
(Johnson noise) decreases (increases) monotonously in mid-frequency range with a. The dip 
in shunt noise in the mid-frequency range is deepest at a « 0.4, and becomes a step-up 
feature for a > 0.6. 
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FIG. 5. Current noise of a two-block hanging model, with varying a = gtes,i/(gtes,i + gtesfi) = 
0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 in linear-log (a) and log-log (b) scales. Increasing a corresponds 
to a shift of all the extra features to higher frequencies (both low- and high-frequency sides), a 
decrease of the phonon noise \I u \tes,b and an increase of the Johnson noise in the mid- frequency 
range, and a non- monotonous behavior of |-/u>|tes,i 

and \I u \sh- The magnitude of \Iu>\tes,i noise 
component and the total noise have a maximum at around a = 0.6, and |iw| s /i a minimum around 
a = 0.4. Other parameters: C = 1.65, Ci/C tes = 10, R = 0.117, p = l,R L = 0.001ft, J = 10/xA, 
g tes ,b = 1 nW/K, T e i = L/[R L + R (l + (3)] = 0.015r tes , T tes = C te s/gtes,b- Dashed line in (a) shows 
the single-block (no C\) result for total noise. 

Finally, we also want to plot the dependence on the loop gain C, shown in Fig. [6j The 
main effect of increasing C is the same as in the simple modeP, which is that it increases 
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the magnitude of the thermal noise components, and reduces the Johnson noise level below 
the frequencies set by the effective thermal time constant. This means that for low values of 
C < 1 the total noise is Johnson noise limited [Figs [6] (a) and (c)], whereas for higher loop 
gains C > 1 the thermal noise dominates [Figs [6] (b) and (d)]. This transition means that 
the relative size of the mid-frequency bump first grows faster with C, and then more slowly 
when Johnson noise becomes irrelevant. Notice that for high £, also the high frequency 
cut-off for the thermal noise components moves to higher frequencies, due to the increase of 
the effective thermal time constant. The shunt noise behaves, again, in a more complex way. 
For low C < 1, the low-frequency part of it is first suppressed like Johnson noise, but for 
high C > 1 the situation reverses, and eventually for the highest values, the low-frequency 
shunt noise is above the high-frequency level. 



2. Intermediate model 



For the intermediate model, equivalent equations to Eqs. 16 can be written and solved 



in frequency domain, with the following simple results for the two TFN noise components: 

I T |2 _ p2 I / A|2 #L,l( T l) 1 



»W 1,6 — 1 1,6 I frp \ , \o 1 i 2 2' 

{gtes,i{Ti) + gi, b ) 2 1 + u 2 t{ 

2 d2 I „ /, \|2 



\Iu>\tes,l — Ptes,l\ S 1(^)1 — H~ 2~T ' ( 2 ^) 

where T\ = Ci/(gt es ,i{Ti) + gi,b), as before for the intermediate model. The Johnson noise 
results (both TES ans shunt) have no direct dependence on the thermal model parameters, 



thus Equations 20 and 21 remain the same, as long as one uses the correct equation for 
Z tes . As Z tes does not qualitatively differ between the intermediate and hanging models, the 
Johnson noise terms also behave qualitatively the same way. 

In Fig. [7] we plot examples of how the noise depends on the most relevant parameters 
for this model. T\ was calculated using the simplifying assumption n = m = 4. The results 
for the total noise as a function of C\ [Fig. [7] (a)] look nearly identical to the results of the 
hanging model, Fig. |4j The only main difference is that the level of low-frequency noise 
is lower because the intermediate block is at a lower temperature than the hanging block, 
thus reducing the TFN noise level. Naturally, the breakdown of the TFN noise into two 
components is completely different in this case: both \I u \i,b and |iu;|tes,i have a low-frequency 
component, with l-Llfe^i developing the bump structure at intermediate frequencies above 
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FIG. 6. (Color online) Current noise of a two-block hanging model, with varying C. In the left 
panels C = 0.33,0.49,0.73, 1.1 in linear-log (a) and log-log (c) scales, whereas in the right panels 
C = 1.65,2.48,3.71,5.57,8.35, in linear-log (b) and log-log (d) scales. Increasing C corresponds 
to an increase of the |/ w |tes,6 an d |-?aj|ies,i TFN noise levels, and an increase of the thermal cut-off 
frequency for C > 1. The low-frequency Johnson noise decreases, and \I u \ s h behaves non- 

monotonously. The minimum of the low- frequency |J w | s ft noise component appears at C ~ 1. Other 
parameters: a = 0.5, Ci/C tes = 10, #0 = O.lfi, = 1, R L = O.OOlfi, I Q = 10//A, g tes ,b = 1 nW/K, 
T d = L/[Rl + Ro(l + /?)] = 0.015r tes , r tes = C tes /gtes,b- Dashed lines in (a) and (b) show a few 
results for the the one-block model for total noise. 

1/ti. Fig [7] (b) shows the dependence on C for high C values. Again, the total noise is 
nearly identical. Interestingly, in the intermediate model the shape of |/ w |te«,i noise actually 
changes with C, by increasing the relative size of the bump-feature. The result of that is 
that the total TFN noise looks nearly the same as in the hanging model, in which the shape 
changes due to the stronger growth of the hanging noise component relative to the phonon 
noise. 
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FIG. 7. (Color online) (a) Current noise of a two-block intermediate model, with vary- 
ing Cx/Ctes = 0.33,0.49,0.73,1.1,1.65,2.48,3.71,5.57,8.35,12.5 in linear-log scale. Increasing 
C\ corresponds to a growth of the bump and a shift of its onset to lower frequencies, and 
a decrease of the cut-off frequency for the |iw|i,6 noise component, (b) Same, with varying 
C = 2.48,3.71,5.57,8.35,12.5. Increasing C corresponds to an increase of the |iw|l,6 and |/ w |tes,i 
noise levels, and an increase of the thermal cut-off frequency for |/ w |tea,i- ( c ) Same, with varying 
a = g t es,i(Ti)/(gte S ,i( T l) + 9i,b) = 0.1, 0.2, 0.3, 0.4, 0.5, and (d) with a = 0.5, 0.6, 0.7, 0.8, 0.9, 0.95. 
Increasing a corresponds to a monotonous increase of the phonon noise level ll^lifc, but a non- 
monotonous behavior of |io;|tes,i [first increase in (c), then decrease in (d)] with maximum bump 
amplitude around a = 0.6 (maximum low- frequency noise around a = 0.4), and a shift of the bump 
to higher frequencies. Johnson noise increases in the mid- frequency range. Other parameters 
in plots (if not varied): C = 1.65, Ci/C tes = 10, a = 0.5, R = 0.1O, /3 = 1, R L = O.OOlfi, 
I = 10/xA, <7i, 6 = 1 nW/K, r el = L/[R L + R (l + /?)] = 0.015r tes , r tes = C tes /g hb . 
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In Figs. |7](c) and (d) we show the dependence on gt es ,i, keeping constant, parametrized 
by the relative strength of the two thermal conductances, a = gtes,i(Ti) / (gtes,i(Ti) + <7i,h)- 
Once again, the dependence is more complicated. The total noise has a maximum at 
a ~ 0.6 — 0.7 at intermediate frequencies (the bump structure), and the bump shifs up in 
frequency with increasing a, and finally vanishes for high enough a, as before for the hanging 
model. The low frequency noise increases because of the increasing effective conductance 
to the bath, in contrast to the hanging case. By looking at the TFN noise components 
separately, once sees that the appearance and disappearance of the bump is entirely due 
to the |iu;|t eS) i component (having a maximum around a = 0.5), as the \I u \i,b part simply 
increases monotonously with a, becoming dominant at low frequencies for high a. The 
Johnson noise again increases in the mid-frequency range with a. For the highest value of a 
here, the total noise already looks like the noise of a single-block model, but with a larger 
heat capacity C tes + C\. 



3. Parallel model 

Finally, for the parallel model [Fig. [I] (c)], due to the one extra thermal link compared to 
the hanging and intermediate models, there will be one more TFN noise component. Solving 
for all the TFN noise components, one gets: 

it 1 2 D 2 I / \|2 9tes,l(Tl) I 

1^1 1,6 = Pl,b\Sl{w 
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{g tes , l {T l )+g 1>b fl + uHr 
2 9tb/(9tes,i(Ti) + 9i,b) 2 + u 2 t? 

1 + L0 2 T 2 

(24) 



with T\ = Ci/(g tes ,i(Ti) + gi,b) (Eq. 13) the same as for the intermediate model, because 
the added gt es ,b does not contact C\. 

In Fig. [8j we show how the different noise components evolve as a function of the 
thermal conductance g^. This means that the starting point is the hanging model with 
gxj, = 0. Again, T\ was calculated using the simplifying assumption n = m = 4. The 
general picture is that the total noise increases as g^b increases, parametrized in the plot as 
b = gi,b/ [gi,b+ gtes,i(Ti)}. In terms of the noise components, the low-frequency shoulder of the 
phonon noise \I u \tes,b decreases, and when b > 0.8, the phonon noise becomes flat up to the 
thermal cut-off frequency. On the other hand, the hanging TFN noise |iuj|tes,i develops an 
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FIG. 8. (Color online) Current noise of the two-block parallel model, with varying b = gi t b/[gib + 
9tes,i( T i)} = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9 in (a) linear-log and (b) log-log scales. Increas- 
ing b corresponds to a decrease of the low-frequency \Iui\tes,b noise, increase of the low-frequency 
l-^ltes l noise and a non-monotonous behavior of the |iu|i,ft noise, with a maximum around b = 0.5. 
Low- frequency Johnson noise increases, and shunt noise has a minimum around b = 0.7. Other pa- 
rameters: £ = 1.65, Ci/Ctes = 10, a = gtes,l/{9tes,\+9tes,b) = 0-5, R = 0.1Q, (3 = 1, R L = O.OOlfi, 
Io = 10//A, 0i,6 = 1 nW/K, T e i = L/[R L + R (l + 0)} = 0.015r tes , r tes — Ctes/gtes,b- 

increasing low-frequency level (as in the intermediate model) quite fast, eventually surpassing 
the low frequency noise contribution from the phonon noise at around b = 0.7 — 0.8. In 
addition, the intermediate frequency bump size also increases with b, and moves to higher 
frequencies. The TFN noise component produced by the direct coupling of C\ to bath, l/Ji^, 
has only a low frequency component, whose cut-off moves up in frequency monotonously 
with b. Also, the strength of the cut-off for |iu>|i,& increases with b, as can be seen from the 
log-log plot, Fig. [8] (b). This is in contrast to the other two TFN noise components. The 
noise level for |iuj|i,6 initially increases, but then reaches a maximum around b = 0.5, after 
which it starts decreasing. For the parameter values chosen here, this last TFN component 
is the smallest, but not insignificant. The low-frequency Johnson noise level increases 
with b, and the shunt noise behaves in a complex manner at low frequencies, where it has 
the highest value at low b, decreases until b = 0.7, and then starts increasing again for 
b > 0.7 [See Fig. [8] (b)]. Finally, in the intermediate bump region the total noise increases 
monotonously with b. However, for low frequencies the total noise initially decreases slightly 
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and has a minimum at b = 0.2 — 0.3 (for the parameter values used), after which the 
low-frequency part also increases. 



III. THREE-BLOCK MODELS 



After the exhaustive discussion of the two-block models, we limit ourselves here to two 
examples of three-block models, which we have already used in analysis of real TES dataP^, 
see Fig. [9j The first model is the analog of the hanging two-block model, where there are 
now two hanging extra heat capacities C\ and Ci- We have named this the 2H model. The 
second model has one intermediate and one hanging block, and we name it the IH model. 
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FIG. 9. Three-block models studied, (a) 2H model, (b) IH model. 



A. Equations to be solved 

The extension of the differential equations defining the problem for two-block models, 
Eq. [T], to the equations for three-block models is straightforward. One more equation is 
generated because of the new heat capacity, and one more power flow term is added, which 
in the case of the 2H model is C(T[ es — Tf) in the equation describing C tes . Then, after 
linearization and transformation to frequency domain, the set of equations become for the 
2H model 

. I Q R (2 + P) T 1 1 1 1 

lUJ-ujes — n ±u) -Lw,tes H J-w,l i -l-u,2 ~r n r w , 

^tes Ti T teS) l Ttes,2 ^tes 
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■ tt T/ [n D /i i OM T ^2H\9tes,l + 9tes,2 + 9tes,b) m /nc\ 

iujLI w = V u - [R l + Rq(1 + p)\ I u j T Uttes , (25) 

Jo 

where all g>j are evaluated at T , r teSj i = C tes /g te s,i, T teSt2 = C tes /g tes ^-, and the definitions of 
the more important time constants t% and the effective loop gain C 2 h are now: 

Ctes 



Tl 



(gtes,l + <7tea,2 + 5 , tes,fc)(l — ^tf) ' 



Tl - 
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Comparing to the two-block equations |4| we see that the first equation has one more 
term due to the added heat capacity C%, the third equation is new, and the last term in the 
equation for I w has a changed coefficient due to the new thermal conductance gt es ,2- 



B. Complex impedance 
1. 2H model 



The complex impedance for the 2H model is calculated from Eqs. 25 the same way as 
for the two-block models, with the result 



Ztes,2H = Ro(l+0) + - C2H r R (2+P) 



d\ 1 &2 1 

1 + itUTj 



1 - C 2 H 1 + lUTi 1 - £ 2 H 1 + iur 2 

(27) 



where we have denoted the relative strengths of the thermal conductances as d\ = gtes,i/(gtes,i+ 
9tes,2 + 9tes,b) and d 2 = 5^,2/ (fiW + 9tes,2 + 9tes,b), and other symbols are defined in Eqs. 



26 By comparing with the result for the hanging two-block model, Eq. [6} we see that 
because of the added C2, a new term appears in the denominator. This term is naturally 
mathematically equivalent with the term for C%, as both C\ and C2 are hanging in this 
model. However, one should bear in mind that the term for C\ is not exactly the same as 
in the hanging two-block model, because the new thermal conductance gt es ,2 also affects it 
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through the pre-factor d\. Also, one sees from Eq. 27 that the effect of C 2 is not additive, 
as the new term is in the denominator. It is therefore not obvious how Z tes 2 H behaves as 
a function of the new thermal parameters C 2 and gt es ,2, and we therefore investigate their 



effect by examples, shown in Fig. 10 



What we observe from Figs. 10 (a) and (b) is that by increasing C 2 , the bulge feature 
caused by C\ first gets smaller, and starts to be pushed to a higher frequency (further right). 
However, when C 2 ~ Ci, the trend turns and the bulge feature at the original place also 
starts increasing again. Importantly, for C 2 < C\ the new term does not necessarily lead 
to a new, additive bulge in the complex plane, but just contributes to the old, existing one 



by changing its shape. As a function of d 2 = gtes?/ (gtea,i + 9tes,2 + 9tes,b) [Fig. [TO] (c)] , 
things look a bit like for parameter a in the two-block case. Interestingly, there is always 
an increased bulge compared to the two-block model on the high frequency side, but also 
on low frequency side for small d 2 < 0.15 for these parameter values. In the limit d 2 — > 1 
Ztes,2H approaches a two-block model with a higher effective TES heat capacity C tes + C 2 . 
The loop gain again has a complex effect, as shown in Fig. [2] (d). In the large loop gain 
limit, the three-block hanging model approaches the simple model, but more slowly than in 
the two-block case. For typical values of £ the effect of the added term is strong, as can be 
seen by comparing the two- and three block results. 



2. IH model 

The starting equations for the IH model look just slightly different from the 2H case, as the 
thermal conductances couple different heat capacities, which can be at different steady state 
temperatures (T for C tes and Ci, and T 2 for C 2 ). The derivation thus follows analogously 
to the 2H model, and we simply quote the end result for the complex impedance. It looks 
mathematically exactly the same as the result for the 2H model, with only some changes in 
parametrization: 



ZtesJH = R (l+(3) + 
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(28) 
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FIG. 10. Complex impedance of the three-block 2H model, with varying parameters. We have 



used Rq = 0.1S1 and f3 = 1 in all plots, and frequencies run typically between 0JTt es = 0.01.. 100 (for 
some curves 10 _4 ..100), with Tt es = Ct es /gtes,b- (a) Z tes ,2H as a function of C 2 , with C2/Ct es = 
0.33, 0.49, 0.73, 1.1, 1.65, 2.48, 3.71, 5.57, 8.35, d/C tes = 10, C = 1.65, a = g tes ,i/(gtes,i + 9tes,b) = 
0.5 and di = gtes,2/(gtes,i+gtes,2+gtes,b) = 0.3. Increasing C2 corresponds to a shrinking of the bulge 
caused by C\, and a growth of a second one at higher frequencies. The two-block model limit (C2 = 
0) is shown as the red dotted line, and the one-block limit (C\ = 0) as the dashed line, (b) Same, 
but with C 2 /C t es = 8.35, 12.5, 18.8, 28.2, 42.3, 63.4, 95.1, 142.7, 214.1, 321.1, 481.7, 722.5. Increasing 
C2 corresponds to a growth of the bulge. The dotted and dashed lines are the same as in (a), (c) 
Z tes ,2H as a function of d 2 (or t/ tes , 2 ), with d 2 = 0.02, 0.06, 0.1, 0.15, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 
C = 1.65, Ci/Ctes = 10, C2/Ctes = 10, and a = 0.5. Increasing d 2 corresponds to the added 
bulge feature moving from the left, low- frequency side to the right, high-frequency side. The 
dotted and dashed lines are the same as in (a), (d) Z tes ,2H as a function of £ = Poa/(gt es ,bTo), 
with C = 0.49,0.73,1.1,1.65,2.48,3.71,5.57,8.35,12.5, Ci/C tes = 10, C 2 /C tes = 10, d 2 = 0.3 and 
a = 0.5. The two-block limits (C 2 = 0) are shown as dotted lines, and some one-block results as 
dashed lines. 26 



where the definitions of the prefactors di, the time constants r« and the effective loop gain 
Cm are now: 



di 
d-2 

Tl 
Tl 
T% 



9tes,l 
9tes,l + gtes,2{To) 

gtes,2(To)gtes,2(T2) 

(SW + gtes,2( T o))(gtes,2(T 2 ) + g 2>b ) 
Ctes 

[gtes,l + gtes,2(T Q )]{\ - £ IH ) ' 

c\ 

9tes,l 

c 2 

gtes,2(T 2 ) + g 2 ,b' 



&IH — t : — - — (rr \~\rj-\ ■ (29) 

gtes,i is always evaluated at T and g 2) b at T 2 , therefore we have omitted the temperature 
dependence from their notation for simplification. g 2 fi{Tb) will only come into play through 



Eq. 18, when calculating the noise amplitude. The true loop gain is again the same as for 
the intermediate two-block model (if one takes into account the notation change C\ — > C 2 ), 
given by £ = P a/(g eff T ), where 

gtes,2 

gtes,2\j-2) + g%b 

Because of the mathematical equivalence between Z teS 2H and Z tes iH, we do not discuss the 



details further here, similar plots to Fig. 10 could be generated 



C. Small-signal responsivity 

Again, derivation of the linear responsivity is straightforward using Eqs. [25] or the equiv- 
alent ones for the IH model. As the model dependent terms can all be lumped into Z tes , both 



three-block models also satisfy Equation 15 In Fig. 11 we plot examples of how \si(u 



behaves for the 2H model, similar plots would also follow for the IH model. The overall 
behavior is the same as with the two-block models, with a partial thermal cut-off moving 
to lower frequencies with increasing C 2 , and a decrease of responsivity above that partial 
thermal cut-off with increasing d 2 = gtes,2/(gtes,i +gtes,2 +9tea,b)- A notable difference to the 



two- block results is that for some values of the parameters C 2 and d 2 [C 2 > C\ in Fig. 11 
(a) and d 2 < 0.6 in Fig. [TT](b)], the responsivity shows a "double knee" structure. 
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FIG. 11. Responsivity of a three-block 2H model, with varying thermal parameters, (a) |s/(w)| as 
a function of C 2 , with C 2 /C tes = 0.33,1.1,2.48,5.57,12.5,28.2,63.4,142.7,321.1,722.5, and d 2 = 
9tes,2/(gtes,i + 9tes,2 + 9tes,b) = 0-5. Increasing C 2 corresponds to a shift of the first partial thermal 
cut-off to lower frequencies and the development of a "double knee" in the intermediate frequency 
range, (b) |s/(w)| as a function of d 2 (or g tes ,2), with d 2 = 0.05,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9, 
C 2 /Ctes = 100. Increasing d 2 corresponds to the decrease of responsivity. We have used i?o = 0.1O, 
P = 1, T el = L/[R L + Rq{1 + P)] = 0.015r tes , Ci/Ctes = 10, a = g tes ,i/(gtes,i + gtesjb) = 0.9 and 
C = 1.65 in all plots (if not varied). Tj es = Ct es / gtes,b- Dashed vertical line shows the simple 
model effective thermal time constant, and dotted curves the corresponding two-block results (limit 
C 2 = 0). 

D. Thermodynamic noise 
1. 2H model 



a. Thermal fluctuation noise The thermal fluctuation noise components for the 2H 
model can be derived analogously to the two-block models, by adding the power fluctua- 
tion terms to equations |25j Straightforward algebra then yields for the three TFN noise 
components 



2^2 



L u\tes,2 



P L.2\si(w) 



l+W 2 ^ 2 ' 

'l+w 2 r 2 ' 
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(31) 



if we assume that the three TFN noise sources are all uncorrelated with each other. The 
power noise amplitudes p are defined as before by Eq. [181 and the time constants as 



T i = Ci/g te s,i, just as for the impedance in Eqs. 26 

b. Johnson noise terms Again, we can derive the TES Johnson noise using the internal 



impedance matrix formulation of section II D 1 Just as in the case of responsivity, the end 
result is that the thermal model dependency is fully accounted for by the TES complex 



impedance, so that Eq. 20 is still valid. The same applies to the external (shunt) Johnson 



noise term, where Eq. 21 can still be used. 



We skip plotting the noise for the 2H model, as Eqs. 31 clearly show that the overall 
picture is analogous to the two-block hanging model. Now, we simply have two independent 



TFN noise bumps, which are multiplied by the responsivity curves shown in Fig. [TTJ Because 
of the monotonously decreasing shape of |s/(w)|, the higher frequency TFN noise bump is 
always suppressed more than the lower frequency bump. 



2. IH model 

For the IH model, the TFN current noise terms are naturally derived similarily, with the 
result 

,\|2 U T l 
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tes,l 


- P 2 1 
~~ r tes,l\ 


% 
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tes,2 


- P 2 1 
~~ r tes,l\ 



' l+U 2 T 2) 

,2 02,6/ '(gtes,2(T 2 ) + g 2yb ) 2 + UJ 2 T 2 



71 l + w 2 r| 

%\h = ^>^\gJXT^+92 fi yi+^ (32) 
where T\ = Ci/g tes ,i and r 2 = C 2 / (gtes, 2(^2) +g 2 ,b), as before for the IH model. The Johnson 



noise for the TES and the shunt follow again from the general formulas Eq. pCfland Eq. 21 



For the IH model, we have plotted a few examples of the noise in Fig. 12 It is naturally 
harder to get as complete picture of the phenomenology as for the two-block models, and 
here we only discuss an example where we keep both gt es ,i and gtes^iTo) constant, and 
assume that the thermal exponents for the thermal conductances connected to C 2 are equal 



n = m = 4. In Fig. 12 (a) we see that an increase of C 2 corresponds to a shift of the 



onset of the total noise bump to lower frequencies, but no added bump amplitude in the 
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FIG. 12. (Color online) Current noise of a three-block IH model, (a) with varying C2/Ct es = 
0.49, 1.1, 2.48, 5.57, 12.5, 28.2, 63.4, 142.7, 321.1, 722.5 in linear-log scale. Increasing C 2 corresponds 
to a shift of the onset of the total noise bump to lower frequencies, an addition of a bump feature 
in \Iui\tes,2 that moves to lower frequencies, a reduction of the |iu;|tes,l noise component, and a 
decrease of the cut-off frequency for the \Iw\2 6 noise component. Mid-frequency Johnson noise also 
increases, (b) Same, with varying C = 2.48,3.71,5.57,8.35,12.5. Increasing C corresponds to an 
increase of the TFN noise levels, and a decrease of the low-frequency Johnson noise, (c) Same, with 
varying c = g2,b/(9tes, 2(^2) + 52.fi) = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9. Increasing c corresponds 
to a monotonous increase of the |iu>|tes,l and |/ w |tes,2 noise components, but a non-monotonous 
behavior of 1 7^ 1 2 6 with maximum amplitude around c = 0.3, and a shift of the cut-off to higher 
frequencies. Johnson noise decreases in the mid- frequency range, (d) same as (c), but shown 
in log-log plot. Other parameters in plots (if not varied): C = 1.65, C\/Ctes = 10, c = 0.5, 
b = gtes,i/(gtes,i +5te,, 2 (Ti))=0.5, R = O.lfi, (3 = 1, R L = 0.001O, Jo = WfiA, <? tes , 2 (Ti) = 1 
nW/K, T e i = L/[R L + R (l + (3)] = 0.015r tes , r tes = C tes / g tes , 2 {T ). 
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total noise for this value of C. An additional bump feature appears in |/ w |tes,2 ; which moves 
to lower frequencies with increasing C 2 . In addition, the amplitude of the hanging |/Jf eSj i 
noise component is reduced, and the cut-off frequency for the | -LI 2,6 noise component moves 
to lower frequencies. Mid-frequency Johnson noise also increases below the frequency of 



the onset of the |/u;|tes,i noise component. The effect of the loop gain is seen in Fig. 12 



(b), with the typical result that all the TFN noise components grow and have a higher 
frequency cut-off with increasing C, and that the low-frequency Johnson noise is suppressed. 
Again, the shape of the total noise changes with increasing C, with | -Z"^ Ues,i becoming more 



dominant over |/ w |t es ,2- In Fig- 12 (c), the effect of increasing g 2 ,ft is studied, by using 
the parameter c = g2,b/ (gtes, 2(^2) + Q2,b)- The result is that the total noise level increases 
and develops a clear bump, whose onset moves to higher frequencies. Both |/ w |t es ,i and 
|L|tes,2 noise components increase monotonously, whereas |/ w |2,6 has a maximum amplitude 
around c = 0.3, with its cut-off moving to higher frequencies with c. Also, the shape of 
the I iu I tea 2 noise changes, with the bump disappearing with high values of c. The Johnson 



noise suppression is strengthened with increasing c. The data in Fig. 12 (c), is also shown 



in log-log scale in Fig. 12 (d), highlighting how the shunt noise is suppressed more with 
increasing c, similar to the Johnson noise. A minimum develops, though, for high values of 
c. 

More examples of the IH model, and particularily in combination with analysing real 
impedance and noise data, are discussed in RefsPSEH. 



IV. CONCLUSIONS 

We have given here a comprehensive discussion of of all the possible variants of two-block 
thermal models and their influence on the complex impedance, responsivity and noise of a 
voltage biased transition edge sensor. The results were derived analytically, and easy-to- 
use formulas were provided for impedance, responsivity and noise. In addition, results for 
two variants of simple three-block thermal models were derived, as well. Example plots were 
generated to show how different parameters of the models affect the observables. In general, a 
more complex thermal circuit reduces the responsivity and increases the noise at intermediate 
frequencies, leading thus to performance degradation in terms of the noise equivalent power. 
The derived theoretical formulas in their current formulation are mainly meant to be used 
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in the analysis of TES experiments, and excellent agreement was already achieved in fitting 
real TES detector dat a! 8 * 9 * 36 -*. However, for detector performance optimization one also needs 
to consider figures of merit such as noise equivalent power, energy resolution, speed and 
electrothermal stability in more detail. Those performance considerations will be possible 
in the future for various kinds of TES detectors with more complex thermal circuits, based 
on the work presented here. 
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